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Abstract 

The dynamics of weak vs. strong first order phase transitions is investigated 
numerically for 2+1 dimensional scalar field models. It is argued that the 
change from a weak to a strong transition is itself a (second order) phase 
transition, with the order parameter being the equilibrium fractional popu- 
lation difference between the two phases at the critical temperature, and the 
control parameter being the coefficient of the cubic coupling in the free-energy 
density. The critical point is identified, and a power law controlling the re- 
laxation dynamics at this point is obtained. Possible applications are briefly 
discussed. 



PACS : 98.80. Cq, 64.60.Cn, 64.60.-i 

* electronic mail: gleiser@peterpan.dartmouth.edu 



Of the many interesting possibilities raised by primordial phase transitions the 
generation of the baryon number of the Universe during the electroweak transition has 
been extensively investigated following the seminal work of Kuzmin, Rubakov, and Sha- 
poshnikov 0. For the purpose of this paper, the important aspect of the electroweak 
phase transition is that it is, in most scenarios proposed so far, a first order phase tran- 
sition. And, at least within the context of the standard model of particle physics, the 
transition is very possibly a weak one; the standard computation for nucleation of criti- 
cal bubbles (see Refs. 0- [|) shows that the thin- wall approximation fails and that the 
bubbles are rather thick ||. 

The possibility that the electroweak transition could be weakly first order has led 
Gleiser and Kolb (GK) to propose a novel mechanism by which such transitions evolve 
jTj. The transition would be characterized by a substantial phase mixing as the critical 
temperature T c (i.e. when the two phases are degenerate) is approached from above, 
followed by domain coarsening below T c . (The slow cooling is provided by the expansion 
of the Universe.) GK modelled the dynamics of this phase mixing by estimating the 
fraction of the volume occupied by sub-critical (correlation volume) thermal fluctuations 
of each phase as a function of the temperature. They neglected the fact that these 
sub-critical fluctuations were unstable and thus found that the system would be equally 
populated by the two phases as it reached T c . The results of Ref. || indicate that GK 
are at least qualitatively correct; there will be a regime in which the transition is weak 
enough that considerable phase mixing occurs even above T c . (It is of course possible, 
although not directly relevant here, that this interesting regime lies beyond the validity of 
the perturbative evaluation of the electroweak effective potential. Presently this question 
does not appear to be resolved 0.) 

Due to the complex nonequilibrium nature of the system, any analytical approach is 
bound to be severely limited. The need for a numerical investigation of this question 
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is clear. This need is even more justified by noting that several of the gross features 
of the electroweak transition may appear in other unrelated physical systems such as 
nematic liquid crystals and certain magnetic materials. Moreover, numerical simulations 
of first-order transitions in the context of field theories (as opposed to discrete Ising 
models ||10|| ) are scarse. Recent work has shown that the effective nucleation barrier is 
accurately predicted by homogeneous nucleation theory in the context of 2+1-dimensional 
classical field theory [|TTJ. These results were obtained for strong transitions, in which 



the nucleation barrier B was large. Nucleation was made possible due to the fairly high 
temperatures used in the simulations. (Recall that the decay time is proportional to 
exp(5/T).) 

In order to study how the weakness of the transition will affect its dynamics, the 
homogeneous part of the (coarse-grained Helmholtz) free-energy density is written as 

W,t) = ^(t 2 -t|)0 2 -|t0 3 + ^ 4 . (i) 

This free-energy density resembles the finite-temperature effective potential used in the 
description of the electroweak transition, where a is determined by the masses of the 
gauge bosons and T 2 is the spinodal instability temperature 0. In the electroweak case 
the order parameter is the magnitude of the Higgs field and the effective potential is 
obtained after integrating out the gauge and fermionic degrees of freedom. Here, we will 
not be concerned with the limits of validity of the perturbative effective potential. The 
goal is to explore the possible dynamics of a transition with free-energy density given by 
Eq. |J and use the results as suggestive of the behavior in the electroweak case. This free- 
energy density is also similar to the de Gennes-Ginsburg-Landau free energy (with the 
elastic constants set to zero) used in the study of the isotropic-nematic transition in liquid 
crystals This transition is known to be weakly first-order; departures from the mean 
field prediction for the behavior of the correlation length were detected as the degeneracy 
temperature is approached from above, signalling the presence of "pre-transitional phe- 



nomena", due to long- wavelength fluctuations observed by light-scattering experiments 
0- 

It proves convenient to introduce dimensionless variables, x = XyfaT 2 , t = 
tyfoTz, A = (j)/y/J2, and 9 = T/T2, so that we can write the Hamiltonian as 



^ 1 = Ud 2 x\ 1 - I V * I 2 +~ (0 2 - l) X 2 - + ±X 4 
9 9 J 2 IVI 2V J 3 4 



(2) 



where a = a/ '(ay/Tz), and A = X/^aT^). (From now on the tildes will be dropped unless 
a new quantity is introduced.) For temperatures above 9\ = (1 — a 2 /4A) _1//2 there is 
only one minimum at X = 0. At 9 — 9\ an inflection point appears at X- m i = a9\/2X. 
Below 9\ the inflection point separates into a maximum and a minimum given by X± = 



ff [l ± yjl - 4A (1 - 1/9 2 ) /a 2 . At the critical temperature 9 C = (1 - 2a 2 /9A)- 1 / 2 the 
two minima, at X = and X + are degenerate. Below 9 C the minimum at X + becomes 
the global minimum and the Ao-phase becomes metastable. Finally, at 9 = 1 the barrier 
between the two phases disappears. 

In order to study numerically the approach to equilibrium at a given temperature 
9, the coupling of the order parameter A with the thermal bath will be modelled by a 
Markovian Langevin equation, 

d 2 X „ „dX dU(X,6) ~, , 

_ = v . J r_,___i_Li + f (Xlt) , ( 3 ) 

where fj = r]/y/aT 2 is the dimensionless viscosity coefficient, and f = £/aT^ 2 is the 
dimensionless stochastic noise with vanishing mean, related to rj by the fluctuation- 
dissipation theorem, {£(x, t)£,(x', t')) = 2rjT5(t — t')5 2 (x — x'). The viscosity coefficient 
was set to unity in all simulations. The lattice spacing was also set to unity in all 
simulations. It turns out that in all cases of interest here the mean-field correlation 
length £~ 2 = U"(X ,9) will be sufficiently larger than unity to justify this choice. In 
future work it would be interesting to see how to generalize the lattice renormalization 



conditions obtained in Ref. O for temperature independent potentials to the situation 



studied here. The Langevin equation was integrated using the fifth-degree Nordsiek-Geer 
algorithm, which allows for fast integration with high numerical accuracy |T^j. The time 
step used was 5t = 0.2, and results were obtained with a square lattice with L = 64. 
(Comparison with L = 40 and L = 128 produced negligible differences for our present 
purposes.) No dependence of the results was found on the time-step, random noise 
generator, and random noise seed. 

The strategy adopted was to study the behavior of the system given by Eq. |fj] at the 
critical temperature when the two minima are degenerate. The reason for this choice of 
temperature is simple. If at 8 C most of the system is found in the X = phase then as the 
temperature drops below 8 C , one expects homogeneous nucleation to work; the system 
is well-localized in its metastable phase. This is what happens when a system is rapidly 
cooled below its critical temperature (rapid quench), so that it finds itself trapped in the 
metastable state. The large amplitude fluctuations which will eventually appear and grow 
are the nucleating bubbles. If, at 9 C one finds a large probability for the system to be in 
the X + -phase, then considerable phase-mixing is occuring and homogeneous nucleation 
should not be accurate in describing the transition. Large amplitude fluctuations are 
present in the system before it is quenched to temperatures below 6 C . For definiteness 
call the two phases the 0-phase and the +- phase. The phase distribution of the system 
can be measured if the idea of fractional area (volume in three dimensions) is introduced. 
As the field evolves according to Eq. |3], one counts how much of the total area of the 
lattice belongs to the 0-phase with X < X_ (i.e. to the left of the maximum), and how 
much belongs to the +- phase with X > X- (i.e. to the right of the maximum). Dividing 
by the total area one obtains the fractional area in each phase, so that fo(t) + f+(t) = 1, 
indepedently of L 2 . 

The system is prepared initially in the 0-phase, /o(0) = 1 and /+(0) = 0. Thus, 
the area-averaged value of the order parameter, (X)(t) = A~ l J XdA is initially zero. 
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The coupling with the thermal bath will induce fluctuations about X = 0. By keeping 
A = 0.1 fixed, the dependence of fo(t), /+(£), and (X)(t) on the value of a can be 
measured. Larger values of a imply stronger transitions. This is clear from the expression 
for 6 C which approaches unity as a — > 0. (In the electroweak case, the same argument 
applies, as what is relevant is the ratio a 2 /A; a is fixed but A increases as the Higgs mass 
increases.) The results are shown in Figure 1 for several values of a between a = 0.3 
and a = 0.4. Each one of these curves is the result of averaging over 200 runs. The 
two important features here are the final equilibrium fraction in each phase and the 
equilibration time-scale. The approach to equilibrium can be fitted at all times to an 
exponential, 

fo(t) = (l - / EQ ) exp [- (t/r EQ r] + / EQ , (4) 

where f is the final equilibrium fraction and teq is the equilibration time-scale. In 
Table 1 the values of teq and a are listed for different values of a. Note that the slot for 
a = 0.36 is empty. For this value of a the approach to equilibrium cannot be fitted at 
all times to an exponential; however, at large times it can be fitted to a power law, 

/(>(*) U=0.36 OC t~ k , (5) 

where k is the critical exponent controlling the approach to equilibrium. A good fit is 
obtained for k = 0.25 (±0.01) as shown in Figure 2. Note that this is not the same as 
the dynamical critical exponent z ~ 2, defined as teq = £,* or - 

The fact that there is a critical slowing down of the dynamics for a ~ 0.36 is indicative 
of the presence of a second order phase transition near a ~ 0.36. Similar behavior has 
been found in liquid crystals in the neighborhood of the isotropic-nematic transition [[15 



and is typical of ferromagnetic materials near the Curie temperature. This transition 
reveals itself in a striking way if we define as an order parameter the equilibrium fractional 
difference AF E q, 



AF EQ = / EQ - . (6) 

In Figure 3 AFeq is plotted as a function of a. Clearly, there is a marked change in 
the behavior of the system around a = a c ~ 0.36. This curve is essentially identical to 
numerical results for the magnetization as a function of temperature in Ising models; the 
rounding is due to finite size effects. (See Fig. 2a in [TO.) For a < ot c the fractional 
area occupied by both phases in equilibrium is practically the same at 0.5. There is 
considerable mixing of the two phases, with the system unable to distinguish between 
them. One may call this phase the symmetric phase with respect to the order parameter 
AF E q. For a > a c there is a clear distinction between the two phases, with the +- 
phase being sharply suppressed. This may be called the broken-symmetric phase. As a 
consequence of this behavior a very clear distinction between a strong and weak transition 
is possible. A strong transition has a > a c so that the system is dominated by the 0- 
phase at 8 C . For a weak transition neither phase clearly dominates and, as argued above, 
the dynamics should be quite different from the usual nucleation mechanism. 

In order to understand the reason for the sharp change of behavior of the system near 
a c , in Figure 4 the equilibrium area-averaged order paramemeter (A)eq and the inflection 



point X in f = ff 1 — J 1 — 3A(1 — I/O 2 ) /a 2 , are shown as a function of a. Also shown is 
the rms amplitude of correlation-size fluctuations A r 2 ms = (X 2 )t — (X)?, where (• ■ -)t is 
the normalized thermal average with probability distribution P[A SC ] = exp [— F[X sc \/9]. 
F[X SC ] is the free energy of a gaussian-shape sub-critical fluctuation. For details see Ref. 
||16|| . It is clear from this figure that the transition from weak to strong occurs as (A)eq 
drops below Xj n f. This result can be interpreted as an effective Ginzburg criterion for the 
weak-to-strong transition. It matches quite well the fact that the critical slowing down 
occurs for a ~ 0.36. This result is in qualitative agreement with the study of Langer 
et al. contrasting the onset of nucleation vs. spinodal decomposition for binary fluid 
and solid solutions |T7|], where it was found that the transition between the two regimes 
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occurs roughly at the spinodal {i.e. at the inflection point). Even though X rms drops 
below X in f for a smaller value of a, being a much less computer intensive quantitity to 
obtain, it should serve as a rough indicator of the weak-to-strong transition. 

The present work raises many questions for future investigation. Apart from in- 
vestigating the 3+1-dimensional case, and obtaining the critical exponent for the order 
parameter (as well as more accurate values for a c , teq and k) using finite-size scaling 
techniques HlOfl , AF EQ oc (a — a c Y ((3 = 1/8 for the d = 2 Ising model, and (3 = 1/2 
for mean field), it should be interesting to test if this behavior could be observed in the 
laboratory. A possible system would be an Ising magnetic film in the absence of an ex- 
ternal field, heated to just above its measured Curie point. For the model studied here, 
this would correspond to a < a c . 
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TABLE 1 



a 




a 




U(0c) 


0.30 


21.0 


0.80 


0.505 


0.495 


0.33 


40.0 


0.80 


0.514 


0.486 


0.35 


75.0 


0.60 


0.525 


0.475 


0.36 






0.580 


0.420 


0.37 


25.0 


0.65 


0.800 


0.200 


0.38 


15.0 


0.80 


0.870 


0.130 


0.40 


5.0 


1.0 


0.937 


0.063 
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